In the diffraction shadow: Norton waves versus surface plasmon-polaritons in the optical region 
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Surface electromagnetic modes supported by metal surfaces have a great potential for uses in miniaturised 
detectors and optical circuits. For many applications these modes are excited locally. In the optical regime. 
Surface Plasmon Polaritons (SPPs) have been thought to dominate the fields at the surface, beyond a transition 
region comprising 3-4 wavelengths from the source. In this work we demonstrate that at sufficiently long 
distances SPPs are not the main contribution to the field. Instead, for all metals, a different type of wave prevails, 
which we term Norton waves for their reminiscence to those found in the radio-wave regime at the surface of 
the Earth. Our results show that Norton Waves are stronger at the surface than SPPs at distances larger than 6-9 
SPP's absorption lengths, the precise value depending on wavelength and metal. Moreover, Norton waves decay 
more slowly than SPPs in the direction normal to the surface. 



The confinement of the electromagnetic field associated 
to Surface Plasmon Polaritons (SPPs), and their intrinsic 
speed, make them very interesting candidates for their use in 
photonics [1 , 2J. Due to this, the study of the electromagnetic 
(EM) fields radiated by localized sources (like defects fS), 
nano-gratings or aperturesf?!) placed on a surface has re- 
ceived a renewed interest in the last decade. This is an old 
problem, which was put at the forefront in the early 1900's 
by its possible relevance to the transmission of radio signals. 
The seminal works of Zenneck|5| and Sommerfeld|6| un- 
veiled the existence of surface waves running along the Earth, 
which can be considered as a lossy dielectric. The interest 
in these works fainted after the realization that radio trans- 
mission does not occur via the exponentially damped sur- 
face modes, but through reflection at the ionosphere. Nev- 
ertheless, Norton subsequently showed that, in the long dis- 
tance limit, radio waves decay algebraically at the surface|7|. 
This result triggered a debate on the range of validity of 
Zenneck-Sommerfeld and Norton waves in the radio regime 
that has propagated to our days (see Ref. fSl for more details 
and an historical account). Recently, advances in nanofab- 
rication have allowed the scaling down of old radio devices 
into the optical regime|9|. Metallic surfaces are specially in- 
teresting because they support SPPs, which are surface EM 
modes strongly confined to the plane. The analysis of the sur- 
face EM fields created by a localized source in a metal sur- 
face has revealed the existence of a near-field region, extend- 
ing for 3-4 wavelengths, where the field presents a complex 
dependence lfTOl [TTl . SPPs have been thought to dominate the 
EM field beyond this region. In this work we show that, irre- 
spectively of the metal considered, the long-distance asymp- 
totic limit of the EM field at metal surface is not the SPP but 
a different type of wave, which we denote as Norton waves 
(NWs) due to their reminiscence to those found in dielectric 
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surfaces. We show the range of validity of SPPs and NWs and 
the distance and field amplitude after which the latter domi- 
nate. 

Although we will show later how the obtained results ap- 
ply to dipole sources, let us concentrate first on the EM 
fields emerging from a subwavelength slit, placed in an op- 
tically thick metal film. The film is back-illuminated by 
normal-incident p-polarized light with wavelength A (i.e. 
the wavevector in vacuum is g = 2n/X). The frequency- 
dependent dielectric constant of the metal is e,„. This system 
has been chosen for analytical simplicity (the full EM field can 
be expressed in terms of the magnetic field along the slit axes, 
Hy{X,Z)) and, also because it is a configuration that has been 
amply studied both theoreticaUy|[l2l [HI [H [15] [H [17 1 and 
experimentallvllTSl [T9l l20l |2T1 . Figure 1 renders a snapshot 
of the radiated Hy{X,Z) (computed with the FDTD method) 
for slit width of A = lOOnm and A = 540«m, for both a Per- 
fect Electrical Conductor (PEC, characterized by \e,„\ = °°) 
and Au|22 |. The choice of metal and wavelength is motivated 
for proof-of -principle purposes on the existence of NWs, but 
we will show later on that our results are applicable to other 
metals and frequency ranges. Our treatment fully takes into 
account the vectorial nature of the EM fields and, therefore, 
goes beyond the scalar approximations considered in other 
works|[T71l20J. 

It is apparent from Fig. 1 that the effect of a finite £„, is 
to strongly modify the radiation pattern close to the metal sur- 
face. Although we will provide expressions for the field every- 
where, our main focus will be to characterize the fields within 
the diffraction shadow which, loosely speaking is the region 
where radiation from a slit in a real metal is strongly reduced 
with respect to the PEC case. 

The Green's dyadic method is more suitable for an analyt- 
ical study of this problem. Within this method, the field radi- 
ated at the point R = {X,Z) by a slit of width A is (see Supple- 
mentary Information SI for the justification of this expression 
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FIG. 1: Snapshot of the magnetic field radiated by a back-illuminated subwavelength slit in an optically thick metal film. In the top panel the 
metal is treated as a Perfect Electrical Conductor, while the metal in the lower panel is Au. The wavelength is A = 540nm. 



and its validation with numerical calculations) 
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H{x,z)~y/e,„-\ I G{x-x',z)E,{x',z = -5)dx' (1) 
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where the "Green's function" G{x,z) is the magnetic field gen- 
erated by a dipolar source with the electric field pointing along 
the x-direction, placed at the metal interface, and 5 is the skin 
depth for the metal. In this expression and throughout the pa- 
per all distances denoted by lower case letters are expressed 
in dimensionless units as x = gX, z = gZ and a ^ gA. Alter- 
natively, given that the fundamental waveguide mode inside 
the slit is constant in the x-direction, G{x,z) can be seen as 
the magnetic field radiated by an infinitesimally thin slit. The 
angular spectrum representation of this function is: 



G{x,z) 



D{q)e 



iqx-\-iq-z 
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where q is the x-component of the wavevector (in units of 
g),D{q) = qzm/[2n{e,nqz + q;„)], q^ = V^^l^ Izm = 



The solution to this integral is not known in closed form. 
Fortunately, there are mathematical methods 1 28 1 for extract- 
ing its long-distance asymptotic expression, Gasymp{x,z)- The 
rigorous calculation for Gasymp{x,z) is provided in the Supple- 
mentary Material S2 and, additionally, a simplified derivation 
will be given later on. But, before going into the mathemat- 
ical details, let us now concentrate on the fields at the metal 



surface and give the result obtained: 

Gasympix,0) = Gsppix,0) + Gnw{x,0) . (3) 
In this expression Gspp{x,Q) is the SPP contribution 

GsPp{x,0)^27:iCpe'''''\ (4) 

where qp — \/£m/i^ + £m) is the SPP momentum, and Cp = 
q^p/ [27r(e,„ ~ 1)] is the residue of D{q) at qp. 



The second term is 



Gnw{x,Q) 



(5) 



As will be shown later, this term is the 2D optical analog in 
metal surfaces of the Norton wave|7| found in the study of the 
radio-wave radiation of point dipoles on lossy dielectric inter- 
faces. Dimensionality accounts for the difference between the 
decay laws: x^^l'^ (2D dipoles) and the x^^ (3D dipoles). 

The validity of Eq. (|3]l and the competition between SPPs 
and NWs is illustrated in Fig. [2] which shows the magnetic 
field at the surface radiated by an infinitesimally thin sub- 
wavelength slit, for Au at two different wavelengths. In each 
case, this figure renders the exact result (computed numeri- 
cally from Eq. ([T])) and the the SPP and NW contributions. 
For the cases considered in this figure, the asymptotic result 
given by Eq. ([3]l is virtually indistinguishable from the exact 
result even for X w 3/im and it is not represented. The field is 
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FIG. 2: The magnetic field at the metal surface radiated by an hori- 
zontal dipole as function of distance. The dependencies are presented 
for Au at two wavelengths: 540 nm (continuous curves) and 567 nm 
(discontinuous curves). The figure shows the exact result (red) and 
the SPP (black) and NW (blue) contributions. 



mainly SPP-like at the shorter distances, while NW dominates 
at sufficiently long distances from the source. Notice that the 
relative phase of the NW and the SPP contributions at the dis- 
tance where their modulus are equal changes with wavelength. 
So, their destructive interference may lead to the cancellation 
of the field (as in Au, at = 540nm at X w I2jj.m) or, if the 
cancellation is not complete, to the appearance of small os- 
cillations in the total field amplitude of the field (as in the 
presented case of Au at A = 567nm). It is worth noticing that 
similar oscillations were found in Scanning Near Field Opti- 
cal Microscope experiments in Au[20|, but their origen was 
unknown. 

Beyond the particular examples presented in Fig. [2j the 
expression given by Eq. (|3]l is a good approximation for 
the field at the surface, at sufficiently long distances from 
the source. In order to quantify this statement, we de- 
fine Xa(j3) as the minimum distance such that |(G(x„,0) — 
Gasymp{xa,0)) /G{xa,Q)\ < j3. The discontinuous lines in 
Fig. [5] show the spectral dependence in the optical and tele- 
com regimes of x„(0.1) for different metals, in units of the 
corresponding SPP absorption length. Given that the NW de- 
cays algebraically and the SPP exponentially with distance, at 
sufficiently large distances the NW is the main contribution to 
the field at the surface, for all metals and all wavelengths. The 
crossover from NW to SPP is represented in Figure. [3] which 
renders the spectral dependence of the distance at which the 
NW contribution is larger than the SPP one, x^w, for differ- 
ent metals. This distance strongly depends upon the dielectric 
permittivity of the metal, being smaller for very lossy metals, 
as Cu and Au in the region of inter-band transition (close to 
A = 500 nm for both metals). 

Undoubtedly, the existence of NWs in metal surfaces has 
passed unnoticed up to now due to their small amplitude. In 
order to characterize how much has the field decayed when 
the NW takes over, we consider the ratio \G{Xnw ,0) / G{X — 
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FIG. 3: Spectral dependence of the crossover between SPP and NW 
at the metal surface, for different metals. The continuous lines render 
xj^iw, defined as the distance at which the amplitude of the NW is 
larger than the SPP one. The discontinuous lines render Xa(O.l) the 
minimum distance at which Gasvmp gives a relative error of 10% with 
respect to the exact result. The inset renders the spectral dependence 
SPP absorption length for the considered metals. 



A/4, 0) I (the distance X — X/4- has been arbitrarily chosen to 
give a representative reference in the near-field). In the op- 
tical regime, when the NW takes over the field has decayed 
by a factor ranging from 10^^ for lossy metals (like Ni) to 
10^^ — lO^'* for Ag. Therefore, the NW is not a good chan- 
nel for sending information along the surface. Nevertheless, 
and given that estimations of the field at the surface far away 
from the source based on the decay of SPPs may be orders of 
magnitude wrong, NWs may have to be taken into account for 
precise analysis or design of experiments. 

In order to show the origin of NWs, and its relation to other 
waves discussed in the literature, as creeping waves and SPPs, 
let us concentrate on the physical interpretation of the field 
radiated by the slit. Additionally, this will lead to a "poor 
man's" (yet correct) derivation of some of the main results. 
It is clear from Eq. (|2]i that a slit excites the whole range of 
diffraction modes (both radiative and evanescent) with an am- 
plitude given by D{q), which can loosely speaking be under- 
stood as the density of EM modes with a given wavevector q at 
the slit position|25 1. The standard treatment of G{x,z) in the 
far-field relies on the observation that, although all modes are 
always present, their contributions cancel out due to destruc- 
tive interference whenever the phase <I> = qr changes rapidly. 
Thus, only the region in ^-space where the phase presents an 
extremum contributes to the far field. For a given point {x,z) 
(or (r, 9) in polar coordinates with 9 defined as the angle from 
the normal to the surface), the extremum occurs at the con- 
dition {q/qz)min = x/z, i.e. ^„„„ = sin0. Expanding the in- 
tegrand around this extremum leads to the "ray-optics" (RO) 
contribution 
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This analytical result reproduces what was observed in 
Fig. 1 : the magnetic field radiated by an infinitesimally thin 
slit in a PEC is isotropic, but the pattern in a real metal is 
strongly modified close to the surface, for angles such that 
COS0 < l/^|e„,|. 

However, right at the surface the derivative of the phase 

= qx never cancels and the saddle point approximation out- 
lined above can not be directly applied. As the integral of 
the product of a smooth and rapidly oscillating function is 
very small, only the parts of the angular spectrum where D{q) 
changes rapidly in the scale of 2n/x will give a net contribu- 
tion to the integral. For very small x all the "density of states" 
contribute. As x increases, the smooth long-^ region of D{q) 
is progressively canceled out in the integral, which is eventu- 
ally dominated by the strong (and rapid) contribution from the 
pole in D{q). The contribution of this pole gives the SPP field. 
Notice that, in a lossy metal, the density of states associated 
to the plasmon pole has a finite width, which causes the expo- 
nential decrease of the SPP amplitude with distance (charac- 
terized by the SPP propagation length Ispp — Im(^p) ^ ' ) . 

The previous argument explains why the field at the sur- 
face is not the SPP for all distances and is expected to have a 
complex dependence with x. Recently, Lalanne and cowork- 
ers have termed "creeping wave" (CW) to the difference be- 
tween the exact field and the approximation given by the SPP 
pole||llJ. The numerical study of the CW has shown that it 
is a damped wave which, along the surface, oscillates with 
the free-space wavevector and decays after a few wavelengths. 
A point to notice is that despite the e" dependence, the CW 
arises from the whole angular spectrum, not only from regions 
close to ^ = 1 . 

However, the SPP pole is not the sharpest feature of D{q): 
the derivative of D{q) diverges at the branch point q, = 0. This 
is illustrated in Fig.|4] which shows that |£>(<?)| has a kink at 
q = I. The contribution to the integral from this kink is ex- 
pected to be small but, as the kink can not be characterized 
by a typical width in q-space, it is not as strongly suppressed 
as the SPP contribution when integrated with an oscillatory 
function. In order to show that the kink originates the NW, it is 
convenient to integrate by parts G(ji:,0). Then, from Eq.(|2| we 
obtain G{x,0) = (i/x) p^A{q)e''>'' dq, with A(^) = qG'{q). 
This representation has the advantage that the kink in D{q) 
transforms into a square root singularity (see inset to Fig. |4|i. 
The contribution close to q^ = can be retrieved by keep- 
ing the singularity but setting q^ — everywhere else, this is, 
by defining A^^ (q) = (e,„/ [2n./e,„^] ) The inset 

to Fig. [4] renders the comparison between A(^) and A'^^lq), 
for a representative case. Of course, A^'^(^) is only a good 
approximation to A(^) close to ^ = 1, so its use for inte- 
gration over the whole angular spectrum could seem unjus- 
tified. However for very large x this is valid, as only the 
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FIG. 4: The modulus of D{q), defined in Eq. (|2]l. The blue and 
red curves are for the Au surface (at A = SOOnm and A = 540nm, 
respectively), whereas the green curve is for the PEC (^ = 0). 
The inset shows the same for the function A{q) defined in the text. 



region close to ^ = 1 contributes. With this, Gnw{x,0) — 

(i/x) r^A^'^{q)e''"^dq = /(£„,/ [2^^;;^])hI^\x)/x. Re- 
calling that this expression is only valid for large x, we sub- 
stitute Hq^\x) for its asymptotic value and obtain the result 
in Eq. (|5]). The motivation for our terminology on this type of 
wave is that the asymptotic term found by Norton, for the case 
of radio waves emitted by a dipole in a dielectric|7 1, also de- 
cays algebraically and originates from the angular spectrum 
close to ^ = 1 . The SPP contribution can also be extracted 
from the previous representation by expanding A^ close to 
q = qp. As the NW and the SPP arise from different parts 
of the angular spectrum, their fields can be directly added up, 
leading to Eq. (|3]l. 

The relevance of NWs with respect to SPPs increases when 
we move away from the surface. Since for z = the NW 
originates from ^-values close to the light-line, its decay with 
distance to the surface is expected to be slower than the ex- 
ponential decay of SPPs. In order to obtain the dependence 
of the Norton wave on both x and z, we have carried out the 
asymptotic analysis of the Green's function given by Eq. ([T]i, 
using the general method described in Ref.|28|. The deriva- 
tion and full asymptotic form can be found in the Supple- 
mentary Material S2. The result up to terms of the order 
^-1/2 vvas already presented in Ref.|16|, where it was shown 
that the long-distance asymptotic expressions are excellent 
approximations even for distances as small as x = 1 (i.e., 
X = X/{27z)). However, the asymptotic expression given by 
Ref |16 | misses some contributions of the order r^^l'^, so it 
can not be used for the problem discussed in this paper. We 
find that the expression for G(x,z) in the far-field (r ^ 1) is 
(see Supplementary Material S2): 

Gasymp{x,z) = Gspp{x,z) + Gro{x,z) + Gnw^X^z), (7) 

where Gspp{x,z) = Gspp{x,0)e"if'-^, with qp, = 
-1/ ^( ! + £„,), Gro{x,z) is given by Eq. and Gnw{x,z), 
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defined as the term that goes as r , is given by 



gir+i7t/4 r ^2 



— D(sin0) cos(0) 
2cos((0-0)/2) 



-3/2 



(8) 



The difference between the exact and asymptotic expres- 
sions AG(x,z) = G(x,z) — Gasymp{x,z) decays as r"^/^. 



To summarize, we have shown that, in the asymptotic limit 
of long distances to the source, the SPPs are not the main 
channel for EM fields at the metal surface. Instead, after a 
few SPP absorption lengths, Norton waves take over This oc- 
curs for any metal and any frequency range. Norton waves 
decay much more slowly than SPPs both along the surface (as 
x^3/2 ^Qj. 2D dipoles) and along the perpendicular direction. 
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FIG. 5: The magnetic field close to the metal surface for Au at 
A = 540 nm. The main figure shows the dependence with z for fixed 
X = 10/im. The exact result G{X,Z) (red curve) and the asymp- 
totic one, Gasympi^ ,Z) (red dashed), are presented, together with the 
contribution from the SPP (black curve), ray optics (green curve) and 
NW (blue line). The inset shows the regions in the X — Z plane where 
the different terms dominate. 
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Figure. |5] presents the comparison 
dependence of the exact G{x,z) and G, 
lOjim. These results show that the asymptotic expression is 
very accurate. Also that, as expected, the NW contribution de- 
cays with distance to the surface much more slowly than the 
SPP one. In both chosen examples the SPP dominates right 
at the surface, but the NW takes over at a finite distance from 
it. However, at sufficiently large z the RO contribution always 
dominates. The inset to Fig. [5] shows, for Au at A = 540nm, 
which of the three "asymptotic" terms (SPP, RO and NW) 
dominates in the X — Z plane. This inset shows that the NW is 
the largest contribution over a "stripe" close to the surface. As 
the region where the NW is dominant satisfies z <C x, the full 
expression for the NW given by Eq. ([8]) can be approximated 
by (for |e„,| > 1) 



\Gnw{x,z)\ ~ \Gnw{x,0)\ |^_^^^^^_y|3 



(9) 



The comparison of this expression with 
tribution allows for an estimation of the 
face at which the crossover between NW 
We obtain znw = IV^ml/ll + \/emzAp 
l/|\/e^|, implies znw ~ |\/e^|, or Znw " 
tice that the algebraic decay of the NW 
this wave arises from the interference of 
ponents. 



that of the RO con- 
distance to the sur- 
and RO occurs, Znw- 
which, if znw/x ^ 
=!A|Ve;;|/(2;r).No- 
with z, reflects that 
its constituent com- 
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APPENDIX A: FIELD IN THE VACUUM HALF-SPACE 
USING THE GREEN'S DYADIC 

Consider a p-polarized electromagnetic wave (with wave- 
length X and wavevector g = 27z/X) incident onto a metal film 
with a subwavelength slit. The metal film is optically thick 
and extends from z = —W (where the EM field impinges) to 
z = (the exit side). The dielectric constant of the metal is e,„. 
The slit has width A and we set the origin of the x-axis at the 
center of the slit. 

According to Lippmann-Shwinger integral equation[26|, 
the electric field at any point at exit side of the film (z > 0) 
is given by the following integral relation 

E(R) =Eo(R)+g^ / t/R'Ae(R')G^(R,R')E(R'), (Al) 



where Eo(r) is the solution without the slit, A£(r) = 1 — e,„ in 
the volume occupied by the slit, V, and zero everywhere else. 

In the case of an optically thick film the field Eo(R) can 
be neglected at the exit side and the dyadic G^(R,R') can 
be approximated by the one corresponding to a single metal- 
vacuum interface. In order to obtain the magnetic field from 
Eq. ( |Al| l, we use the Maxwell equation H = {—i/g)^R x E 
and arrive at 



H(r)= / dr'Ae(r')G^(r,r')E(r'), 



(A2) 



where we have passed to dimensionless distances r = gR. 
The dyadic G^ connects the magnetic field outside of the slit 
with the electric field inside the slit. For the considered two- 
dimensional geometry, where only p-polarized waves are in- 
volved, the magnetic field H points along the y-direction. As- 
suming that the electric field inside the slit mainly points along 
the x-direction only the yz element of the dyadic G^ needs be 
computed. We denote this element by G^(jc,z;x',z'). Follow- 
ing Ref . |27j we find that 



G%{x,z;x\z') = / dq 

In J £mqz + qzm 



(A3) 

where q is the dimensionless .jc-component of the wavevector, 
q = kx/g, qz = \/l-q^ and q, = ^e^-q^. 
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The integrant contains the exponential factor e^"?-"'^ , which 
decays in the distance of a skin depth 5=1 /lm{q^„,), which is 
of the order of a few tens of nm in the optical regime. There- 
fore, the integration limits in z' can be extended to [— oo,0]. 
Moreover, the variation of the Green's dyadic is much faster 
than that of the electric field inside the sUt, so the electric field 
inside the slit can be approximated by its value at the distance 
z = —5 (this is obtained as the average distance to the surface, 
weighted by the exponential decay of the field). An additional 
advantage of using the field at a short distance inside the slit is 
that the numerical problems related to the treatment of corners 
are eliminated. 

The integration over z' can be performed in the following 
way 



slit 



We then obtain 



Izni 



(A4) 



a/2 

-a/2 



with 



1 f 1 

Gsii,ix,z) = — / dq- ■ 

2n J e,„qz + qzn 



iqx+iq~z 



(A5) 
(A6) 



Since we are interested at the field close to the surface and 
away from the source (which, as we will show, arises from the 
region in the angular spectrum close to ^ = 1) and for optical 
frequencies (where £,„ is large), this result can be related to the 
field radiated by a dipolar source G{x,y) at the metal surface: 

l-a/2 

Hy{x,z) ~ \/ e,„ — 1 / dx'G{x — x',z)E^i{x',z = —5), 

J-a/2 

(A7) 

where 



G{x,z) = 



27t 



dq ^ 



(A8) 



The relation between Eq. jAS) and Eq. ( |A7| i can be easily 
seen by noticing that qyj„ « y/ £,„ — 1 in the region of interest. 

While Eq. ( |A7[ ) involves an additional approximation, we 
have preferred to work with the Green's dyadic for a dipole 
source (rather than with Gsih) due to its wider applicability to 
other problems. In any case, the differences for the slit case 
are minimal and the methods described in this paper could be 
straightforwardly applied to Gsit- 

In order to validate the expression ( |A7[ ), we first performed 
full-vectorial computations using the Finite Element Method 
(FEM) of the magnetic field emerging from slits with the 
thicknesses 50 nm and 600 nm. The wavelength was cho- 
sen to be 576 nm and the metal of the film is gold. Then we 
computed the integral given by Eq. ( |A7| i extracting the elec- 
tric field inside the slit at z = —5 from the FEM calculations. 
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FIG. 6: Comparison between the FEM and the analytic computation 
for the magnetic field on the gold-vacuum interface emerging from a 
slit. The slit widths considered are 50 nm and 600 nm, and the wave- 
length is 567 nm. The dashed lines are for the FEM calculations, 
while the continuous lines correspond to the approximate result us- 
ing Eq. jATji. 



The function G{x,z) was substituted by its asymptotic value, 
see the next section S2. The comparison is shown in Fig.|6] 

These results clearly show that Eq. ( |A7| l is very accurate for 
subwavelegth slits (A << A) and even provide a good approx- 
imation for A ^ X. 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF THE FIELD: 
STEEPEST DESCENT METHOD 

Consider the Green's function for a dipole placed at the 
metal surface given by Eq. ( |A8| l. In this section we sketch 
the asymptotic analysis performed, which has been done fol- 
lowing the general method described in lESl for treating Som- 
merfeld integrals. 

In this method, the integrant is first prolonged into the com- 
plex ^-plane. Subsequently, the following changes of variable 

are performed q = sm<j) and s = \/2e' ? sin ^ . In polar 

coordinates (x = rsin0 and y = rcos 9) the integral takes the 
following form 



G = e"' / ds'P{s)e- 



(Bl) 



with 



<I>(^): 



1 Vle-"'/'^ cos[<j){s)]Je,„-sm^[(j){s 



2n 



cos 



£,„ COS[0(.?)]- 



(B2) 

where the contour C in the i-plane corresponds to the real axis 
in the initial plane q. 

Then the integrant is separated into singular and non- 
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FIG. 7: The poles q = ±qp and the branch cuts Im{q^) = with 
infinitesimally small absorption of the vacuum in the complex plane 
q. The poles correspond to the dots. For x > only the pole q = qp 
contributes to the result. 
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FIG. 8: The contours of the integration in the complex plane s for 
different angles 6. The positions of the poles Sp are mark the dots. 
When the pole is crossed during the transformation of the initial path 
to the steepest-decent one, the residue must be taken into account. 



singular parts 



with 



C„ 



<^{s){s~Sp)~Cp 



where Cp is the residue given by Cp 



(B3) 



(B4) 



the position of the pole in the complex plane s is Sp — 
^2e"^/'*sin(^i^), where (j)p = arccos ^ y^qr^ ^ Then, after 
deforming the integral path to the steepest descent one (real 
axis in the plane s), the singular part yields the complemen- 
tary error function inCpe''^'^^i'^ eifc{~is py/r) with the argu- 
ment being the square root of the "numerical distance" in- 
troduced by Sommerfeld. The non-singular part of the inte- 
gral can be expanded in Taylor series close to the saddle point 



s ~0 providing the infinite sum of the integrals of the Gaus- 
sian type. The final result reads 



G = inCpe'^'ip erf c{~ is p^fr) + e"' ^ 



nGeven «!r 2 



ds" 



\s=0- 



(B5) 

The function <I>o is composed of two parts: <t> and Cp/{sp — s). 
Note that the part of the sum in ( |B5[ ) coming from the term 
Cp/ {sp — s) coincides, up to a sign, with the asymptotic expan- 
sion of the complementary error function for large arguments 
without the first term (while the last term coincides with the 
residue contribution into the integral). This expansion reads 



erfc(-/sp Vr) = 20(0 - Op) + i- 



nGeven [Spr] 2 



(B6) 



where © is a Heaviside step function, and Qp is the 
angle defining the diffraction shadow, Qp = Re(0p) — 
arccos (l/cosh[lm(0p)]). This is the critical angle such that 
for Q > Qp^e have Im(sp) < and the initial transformation 
of the integration path into the steepest-decent one leads to 
the crossing of the pole, so that the residue must be taken into 
account, see Fig. |7] An example of the critical angle limiting 
the diffraction shadow is represented in Fig. 1 of the article 
by dashed lines. 

Retaining the residue contribution, Gspp, and the first two 
terms G^o (proportional to r^'/^) and Gj^y^ (proportional to 
r^^/^) in the sum coming from <I> we obtain an expression 
that is correct in the far-field, up to terms AG = 0(r^^l^). 
The result can be rewritten as 



G — Gspp + Gro + Gnw + AG. 
The plasmonic term is 

GsPP = 2niCpe'"i<', 
the "ray optics" term is 

^ir-iK/A coseVem-sin^e 



Gro — 



^2nr e,„ cos Q + V £m - sin^ Q ' 
the Norton wave is given by 



(B7) 



(B8) 



(B9) 



1 e"-'T 



I cos (p\/ e — sin ^ 



cos( 2 ) ecos(^ + Y e — sin^ 

(BIO) 



and and the rest is 



AG = iTTCpe'''-'--'!'^ [erfc(-«px/7) - 2]- 

r(i±^) / d"^is) d" Cp 



MGeven n\r 2 



ds" ds" s — s 

C7„ = 



P/ s=0 

1, n>4, 
0, n < 4. 
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